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PRECISE  INTERFEROMETRIC  PHASE  DETERMINATION 


I .  INTRODUCTION 

This  report  is  primarily  concerned  with  precise  phase-shift 
determinations  where  the  phase  changes  rapidly  (steep-gradient 
region)  or  where  the  phase  changes  are  small  (outer,  low-density 
region).  The  medium  is  assumed  to  have  cylindrical  symmetry. 

Two  things  are  required:  A  precise  knowledge  of  the  initial 
phase  (from  a  pre-shot  interferogram)  and  a  precise  determination 
of  the  final  phase  by  interpolating  the  shifted  fringes  of  the 
data  interferogram.  Abel  inversion  of  the  phase-shift 
information  may  reveal  important  low-density  structures  which 
would  be  missed  without  using  these  techniques.  One  should  note, 
in  the  Abel  inversion  procedure,  that  the  density  determination 
at  a  particular  cylindrical  radius  does  not  depend  on  what 
happens  at  smaller  radii.  Non-uniform  fringe  spacing  may  be 
important  in  the  steep-gradient  regions  and  is  included  in  the 
interpolation  procedure. 


Manuscript  approved  December  1,  1987. 
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The  purpose  of  this  report  is  to  provide  a  general  background 
and  to  set  forth  the  method  in  a  step-by-step  procedure  that  can 
be  put  to  practical  use.  The  book-keeping  details  turned  out  to 
be  rather  involved  and  tricky.  Thus,  it  seemed  reasonable  to 
include  a  detailed  procedure  for  obtaining  the  necessary  input 
from  the  pair  of  pre-shot  and  data  interf erograms  as  well  as  a 
listing  of  the  computer  program.  The  procedure  is  applied  to  a 
sample  interferogram  pair  (Fig.  (1))  obtained  from  a  laser- 
produced  plasma.  The  procedure  includes  preparing  the 
interferogram  pair  for  data  extraction,  tabulation  of  the  data, 
and  input  to  the  computer  program.  The  computed  results  (radius 
and  phase  shift)  are  given  for  the  sample  interferogram  pair.  l 
have  attempted  to  make  the  program  readable  by  the  liberal  use  of 
comments  to  explain  and  de-limit. 

The  example  illustrated  in  Fig.  1  is  an  experimentally- 
obtained  interferogram  pair  from  a  plasma  containing  a  laser- 
produced  blast,  wave  (I~9  to  16).  The  method  to  be  described 
provides  useful  phase-shift  information  for  the  region  of  the 
blast  wave  and  the  outer-most,  low  density  region  (l>16). 

However,  the  central  region  (I<7)  shows  such  irregular  phase 
variation  that  the  method  offers  no  advantage  in  that  region. 
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Symmetry 
Axis  (R=0) 
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BASIC  METHOD 


II . 

We  wish  to  determine  the  radial  density  profile  at  a  given 
axial  postion  Z,  which  is  represented  by  the  line  Z  -  constant 
(called  the  "cut")  in  the  R-Z  plane  of  the  interferogram  pair. 

See  Fig.  1.  The  final  (shifted)  fringes  in  the  region  of  the  cut 
are  first  divided  into  a  sufficiently  large  number  of  radial 
"intervals"  that  the  fringes  are  linear  in  each  interval.  Next, 
the  region  within  the  interval  and  between  the  two  fringes 
bounding  the  cut  is  divided  into  10  regions  by  introducing 
decimal  "fringelets"  which  may  be  spaced  non-uni formly  to 
accomodate  information  in  differing  adjacent  fringe  spacings. 

The  fringelets  are  illustrated  in  Fig.  2  with  dashed  lines  but, 
in  practice,  are  only  treated  numerically  with  the  computer 
program.  The  non-uniform  interpolation  or  fringelet  spacing  is 
described  in  the  appendix.  Each  fringelet  carries  a  Itnown  final 
phase.  The  intersections  (radial  position)  of  the  fringelets 
with  the  cut,  thus,  have  a  known  final  phase.  If  no  fringelet 
intersects  the  cut  within  the  interval,  the  mid-radius  of  the 
interval  is  assigned  the  arithmetic  average  of  the  phases  of  the 
two  fringelets  which  "straddle"  the  cut.  Each  intersection  or 
straddle,  which  assigns  a  known  final  phase  to  a  given  radius 
within  the  interval,  is  called  a  "hit".  The  phase  shift  at  a  hit 
is  obtained  by  subtracting  the  initial  phase  at  the  this  radius. 
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The  initial  phase  is  determined  by  the  phases  of  the  bounding 
initial  fringes  and  by  interpolation  using  the  fractional 
distance  (at  the  hit  radius)  of  the  cut  from  the  lower  fringe. 

The  implementation  of  the  calculations  involves  a  bit  of 
bookkeeping  and  indexing,  some  of  which  is  discussed  here.  Phase 
information  is  recorded  by  sequentially  numbering  the  final  and 
initial  fringes.  They  must,  of  course,  merge  to  the  same  values 
in  the  outer  most,  lowest  density  region  of  the  medium.  The 
initial  fringes  are  usually  close  to  straight  lines  but  are 
divided  in  the  region  of  the  cut  into  a  small  number  of  radial 
"spans"  in  which  the  fringes  are  accurately  linear.  See  Fig.  3. 
We  are  primarily  interested  in  the  pairs  of  initial  and  final 
fringes  which  bound  the  cut  at  a  given  radius.  However,  due  to 
fringes  crossing  the  cut,  these  pairs  change  with  radius.  Each 
interval  II  determines  the  span,  with  index  K  (II),  and  the 
crossing  region  (for  final  fringes),  with  index  L  (II)  of  that 
interval.  In  turn,  K(I1)  and  L(I1)  determine,  respectively,  the 
bounding  initial  and  final  fringe  pairs  for  that  interval.  The 
phase  of  the  lower  bounding  fringe  in  an  interval  is  determined 
by  the  fringe  number;  Ml(K(Il))  for  initial  fringes  and  N(L(I1)) 
for  final  fringes.  Each  hit  (intersection  or  straddle  of  a 
shifted  fringelet  with  the  cut),  with  index  I,  is  labeled  with 
the  interval  number  Il(I)  and  fringelet  number  Jl(I).  See  Fig. 

2.  A  variable  E(l)  indicates  whether  a  hit  is  an  intersection 
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(E(I)-0)  or  is  a  straddle  (E(l)-.5),  The  interpolated  shifted 
fringe  number  at  the  hit  is  thus  N(L(I1(I)))  +  ( Jl  ( I ) -E(  I )  )/lj3'. 
The  interpolated  initial  fringe  number  at  the  hit  depends  on  the 
fraction  F(I)  into  which  the  cut  divides  the  fringe  spacing  at 
the  hit  radius.  See  Fig.  3.  Here,  F(I )— Wl( I )/{W2( I )-(Wl( I ) ) . 
Thus,  the  interpolated  initial  fringe  number  is 
Ml ( K( ll ( I ) ) ) +F( I ) .  The  phase  shift  (in  units  of  pi)  is  thus 
given  by  twice  the  difference  in  the  shifted  and  initial  fringe 
number  at  the  hit. 

P(I)=2*(N(L(I1(I)  )  )  +  ( J1(I)-E(I))/10-(M1(K(I1(I)  )  )+F(I)  ) ) 

The  final  result  is  the  hit  radius  R(I)  and  phase  shift  P(I), 
sorted  for  increasing  radius. 
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III.  PROCEDURE  AND  EXAMPLE 


A  sample  interf erogram  pair,  prepared  for  extracting  data,  is 
shown  in  Fig.l.  This  is  the  superposition  of  two  separate 
interf erograms :  The  inital  or  pre-shot  interferogram  (shown 
dashed)  and  the  final,  shifted  or  data  interferogram.  The 
invidual  interferograms  were  projected  and  the  fringes  carefully 
aligned  to  each  other  for  reference  positions  and  merging  to  the 
same  outer  fringes.  This  happened  to  be  the  interferogram  pair 
fror  which  the  procedure  was  developed  and  was  not  the  optimum 
pair  of  interferograms  for  accurate  application  of  the  method. 

It  may  be  helpful  to  refer  to  Figs.  2  and  3,  showing  separately 
information  on  the  initial  and  final  interferograms. 

We  first  draw  in  the  cut  at  the  desired  Z-position,  taking 
into  consideration  the  requirement  of  cylindrical  symmetry.  We 
next  divide  the  initial  fringes  into  a  small  number  (4  ,  in  the 
example)  of  radial  regions  or  spans  such  that:  (1)  Initial  fringe 
crossings  only  occur  at  jpan  boundaries  and  (2)  the  fringes  are 
accurately  linear  within  a  span.  The  symmetry  axis  (R=0)  and 
outer  radius  (of  fringe  shifting)  are  taken  as  span  boundaries. 
Spans  are  designated  (K)  by  their  left,  or  inner  boundary.  Span 
boundaries  (vertical  lines)  are  extended  into  the  region  above 
the  cut  and  the  following  symbols  used  to  aid  recognition: 

(for  span  boundary),  L  (for  initial  fringe  crossing  the  cut). 
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Each  span  is  labeled  with  the  span  index  (K)  and  lower  initial 
fringe  number  (Ml)  of  fringes  bounding  the  cut.  It  is  also 
helpful  to  note  the  span  boundary  radii  at  each  span  boundary. 

Information  about  the  shifted  fringes  is  indicated  in  the 
region  well  below  the  cut.  First,  locate  all  shifted  fringe 
crossings,  pass  a  vertical  line  through  each  crossing  and 
indicate  with  a  slanted  dash  (^)  that  it  is  a  crossing.  Label 
each  crossing  region  (below  the  cut)  with  the  crossing  index  L 
and  lower  shifted  fringe  number  N  e.g.,  (  ^  n-4  't  ^  '  Next, 

draw  in  (short  vertical  lines)  interval  boundaries  such  that: 

(1)  all  fringe  crossings  (initial  or  shifted)  and  span  boundaries 
are  interval  boundaries,  and  (2)  the  final  (shifted)  fringes  are 
linear  in  each  interval.  Intervals  are  designated  (I)  by  the 
left  or  inner  boundary.  The  symmetry  axis  (1=1,  where  R=0)  and 
outer  radius  (I  =  Ij?+l,  where  10  is  the  number  of  intervals)  are 
also  interval  boundaries.  The  interval  number  of  the  shifted 
fringe  crossings  are  indicated  below  the  crossing  (e.g.,  1=7). 

We  are  now  in  a  position  to  organize  the  information  in 
tabular  form,  on  data  sheets,  so  that  it  can  be  input  to  the 
program.  At  the  top  of  the  first  sheet,  note  the  identifying 
(e.g.,  shot)  number,  the  magnification  (needed  in  Abel  inversion) 
and  the  number  of  intervals  (I^),  spans  (K0)  and  shifted  fringe 
crossings  { N/i )  .  Since  lower  fringe  numbers  (used  in 
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calculating  phase)  depend  directly  on  the  span  (initial)  and  the 
crossing  region  (shifted)  while  information  is  input  for  each 
interval,  it  is  useful  to  relate  each  interval  index  to  a  span 
and  crossing  index.  This  cross-indexing  is  calculated  by  inputs 
from  Table  I.  Next  we  record  (in  Table  II)  the  lower  fringe 
numbers . 

The  next  step  is  to  make  detailed  measurements  of  radial  and 
axial  distances  and  to  organize  this  data  into  tabular  form.  In 
the  shifted  fringes,  it  is  necessary  to  measure  and  record  (in 
Tables  III  and  IV),  for  each  interval,  the  interval  index  I,  the 
radial  position  R1 ( I )  of  the  left  boundary,  the  axial  position 
(negative)  Zl(I)  of  the  lower  fringe  from  the  cut,  the  fringe 
slope  S0( I )=(dZ/dR)/ABS(dZ/dR) ,  and  the  lower,  central,  and  upper 
fringe-pair  spacings  Dl(I),  D2(I),  D3(I).  The  distances  are  all 
measured  in  centimeters.  The  fringe  slope  is  ^  for  the  cut  and 
+1  or  -1  according  as  Z  increases  or  decreases  with  R  for  a 
fringe.  CAUTION ;  At  a  fringe  crossing,  we  can  have  either  a  J=0 
intersection  (where  Zl(I)*(^,  since  St^(I)=-l)  or  a  J=10 
intersection  (where  Zl ( I ) “-D2 ( 1 ) ,  since  S0(I)=1). 

The  shifted  fringe  input  is  first  recorded,  in  Table  III  for 
the  outermost  radius  (where  1=10+1).  The  shifted  fringe  input  at 
all  left  interval  boundaries  (1=1  to  10)  is  recorded  in  Table  IV. 


II 


The  input  data  for  the  initial  fringes  are  recorded  at  each 
span.  This  includes  the  radius  R2{K)  of  the  (left-hand)  span 
boundary,  the  axial  positions  (negative)  of  the  left-side  Z2(K) 
and  right-side  (Z4(K)  of  the  lower  fringes  and  the  axial 
positions  (positive)  of  the  left-side  Z3(K)  and  the  right-side 
Z5(K)  of  the  upper  fringes.  The  values  are  recorded  (Table  V) 
for  the  left  boundaries  of  the  individual  spans. 

The  detailed  shi f ted-f ringe  information  for  each  interval 
(Table  IV)  is  entered,  before  the  program  is  run,  into  DATA 
statements  near  the  end  of  the  program.  The  data  for  cross¬ 
indexing  (Table  I),  lower  fringe  identification  (Table  II),  outer 
radius  information  for  shifted  fringes  (Table  III),  and  the 
intial  fringe  information  (Table  V)  ate  entered  during  the 
running  of  the  program,  as  requested  by  input  statements. 


Interferometry  Data  Sheet  -  1 
(Enter  as  requested  by  program,  during  run) 


Shot  No.  14760  Mag.  M-19 

Ipr.21  KJ?-4  N0r=4 


Table  I.  Cross-Indexing;  The  intervals  (left  interval 
boundaries)  12  which  occur  at  interior  fringe  crossings  and  the 
intervals  (left  boundaries)  13  which  occur  at  interior  span 
boundaries  are: 


Crossing  Index 
L-1 

Interval  Index 
I=I2(L-1) 

Span  Index 
K-1 

Interval  Index 
I=I3(K-1 ) 

1 

7 

1 

3 

2 

10 

2 

7 

3 

15 

3 

18 

Table  II.  The  lower  (shifted/initial 
(region/span)  index  ( L/K ) . 

)  fringe  numbers  vs. 

Index 

Fringe  # 

Index 

Fringe 

L 

N(L) 

K 

M1(K) 

1 

3 

1 

1 

2 

4 

2 

2 

3 

5 

3 

3 

4 

4 

4 

4 

Table  III. 

Shifted  fringe  input  for 

the  outermost  radius 

(I!?+1) 

Rl-15.3  Zl-0  S0=-1 

Dl-  .30  D2-.30  D3=.30 
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Interferometry  Data  Sheet  -  2 
(Enter  in  Data  Statements;  1  interval  per  statement) 

Table  IV.  At  the  left  (inner)  boundary  of  each  interval,  the 
index  I,  radius  Rl(I),  axial  position  of  lower  fringe  Zl(I), 
fringe  slope  S0(l),  and  lower,  central  and  upper  fringe-pair 
spacings  Dl ( I ) ,  D2(I),  D3(I)  are: 


I 

R1(I) 

21(1) 

S^(I) 

Did) 

D2(I ) 

D3(I) 

1 

0 

-.19 

-1 

.40 

.40 

.35 

2 

0.87 

-.22 

-1 

.38 

.38 

.38 

3 

1.78 

-.30 

-1 

.40 

.40 

.40 

4 

2.52 

-.26 

-1 

.38 

.  38 

.  32 

5 

3.69 

-.23 

-1 

.36 

.38 

.34 

6 

5.02 

-.25 

-1 

.36 

.35 

.  30 

7 

6.50 

0 

-1 

.37 

.32 

.32 

8 

7.10 

-.08 

-1 

.38 

.  30 

.32 

9 

8.40 

-.14 

-1 

.  38 

.  30 

.28 

10 

9.02 

0 

-1 

.28 

.30 

.  30 

11 

9.28 

-.02 

-1 

.30 

.32 

.33 

12 

9.45 

-.07 

-1 

.31 

.33 

.34 

13 

9.67 

-.13 

0 

.27 

.33 

.34 

14 

9.90 

-.12 

1 

.24 

.30 

.34 

15 

10.18 

-.25 

1 

.20 

.25 

.25 

16 

10.37 

-.14 

1 

.20 

.24 

.33 

17 

10.69 

-.10 

1 

.29 

.22 

.32 

18 

10.92 

-.10 

0 

.30 

.22 

.33 

19 

11.51 

-.13 

-1 

.32 

.27 

.29 

20 

12.52 

-.17 

-1 

.  30 

.28 

.29 

21 

13.58 

-.20 

-1 

.29 

.29 

.  30 
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Interferometry  Data  Sheet  -  3 
(Enter  or  requested  by  Program,  during  run) 


Table  V.  Initial  Fringe  Information.  The  Radii  and  Axial 
Position  at  the  Span  Boundaries  are: 


Span  Index 

K 

Radius 

R2(K) 

Lower 

Z2(K) 

Fringe 

Z4(K) 

Upper 

23(K) 

Fringe 

Z5{K) 

1 

0 

-.18 

-.30 

.12 

0 

2 

1.78 

0 

-.31 

.31 

0 

3 

6.50 

0 

-.30 

.30 

0 

4 

10.92 

0 

-.29 

.29 

0 
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.  Computer  program  with  data  input  for  the  interferogram 
pair  shown  in  Fig.  1. 

13  r’RUJT  ■r4[F0I.  ft  PROGRftri  USING  NCN-UNIF3Rf1LY-SPftCE3  INTERFERENCE  FRINGE  DECI- 
MftL-  INTERPGI-ftTrON.  OECIfIflL  ( INTERPOLftTEO )  FRINGES  ftfic  CALCEa  FSINSELET3.' 

20  PRINT  -THE  RADIAL  UARIATION  OF  PHASE  IS  QETERNINED  AT  A  FIXED  AXIAL  POSITION 
OR  LINE,  CALLED  THE  CUT.* 

23  INPUT  ‘SHOT  OR  IDENTIFYING  NUN8ER  AND  MAGNIFICATION' i  ID  .fIG 

43  PRINT  -THE  SHIFTED  FRINGES  ARE  MEASURED  FOR  MANY  RACIAL  INTERL'ALS  WITH  IN 
DEX  I- 

50  PRINT  -THE  INITIAL  FRINGES  ARE  TAKEN  AS  LINEAR  FOR  A  FEU  RADIAL  SPANS 
WITH  INDEX  K! I )' 

SO  PRINT  -ALL  SHIFTED  AMO  INITIAL  FRINGE  CROSSINGS  ARE  TAKEN  AT  INTERUAL  BOUN¬ 
DARIES. ALL  UNSHIFTED  FRINGE  CROSSINGS  ARE  AT  SPAN  BOUNDARIES.' 

70  INPUT  -THE  NUMEER  10  OF  RAOI.AL  INTERUALS  IS' i  10 

30  INPUT  -THE  NUMBER  K3  OF  RADIAL  SPANS  IS'iKO 

30  INPUT  -THE  NUMBER  NO  OF  SHIFTED  FRINGE  CROSSINGS  IS'iNO 

100  DIM  A(  SO  )  .3:  33  1  ,C:  20  )  .Dl  (301  .32(  30  )  .D3(  30  )  .£•.  SO  )  .G^  30  )  .H(  SO  ;  .  1 1  ( E3  !  .  13;  1 3  )  .  I 
3(10  )  ,JU53  1  .KiSO  )  ,L'  30  ;  ,MU  10  )  ,N!  13  1  ,R1  (301  .RZC  1 0  )  ,S0(  30  )  ’. 2 H  33  ) 

110  DIM  :3( 10  )  ,23(  10  ) ,24; 10  )  ,2S( 13) 

133  DIM  U(  S3 , 1  i  )  ,X'. 53 , 1  1  )  .Xl  (S3 , 1  1  ( .Yl 53 . 1 1  )  .2; S3 . 1  I  )  .23;  £3 . 11  ) 

133  ■ - L  — - - - - 1 - 3 - 3 - — - - - - 

ua  - Shifted  fringe  calculaticns - 

ISO  - - - 

ISO  ■ - - ---INGEXING  TO  INTERUAL  AND  MUMEERInG-  — - - 

173  PRINT  •intervals  ARE  LABELED  BY  THEIR  INNER  OR  LE-T  BOUNCARY' 

1S3  PRINT  'THE  SHIFTED  FRINGE  CROSSINGS  ARE  LGCATED  AT  INTERVAL  BOUNDARIES  lO+l 
AND  AT 

153  ■ . I  >  CROSSING  INCE.<  L-:  FOR  FOLLDU'NG  LOOP - 

233  FOR  I=:  TO  N3-1 

310  I'iFuT  •:;■  --1  >•  ii:- 1 ) 

233  NeXT  I 
333  L«3 

243  ■— .  I>INTESV'Al  index  FOR  FCLLGUING  LOOP . - 

2S3  Ul  -1 

263  FDR  '=■!  TO  ICINO-I  ) 

270  IF  i;:2(UI  1  then  290 
280  Ul-Ul ♦! 

290  L< I )-UI 
330  NEXT  I 

313  FOR  I-i:(N3-l  )  TO  IS 
320  L< I  )-N0 
330  NEXT  I 

340  print  -the  shifted  (LCUER)  fringe  NUMSER9  APE' 

3E0  FOR  L-l  TO  N3 
563  INPUT  'Nl  L  )-■  iN',  L  ) 

370  NEXT  L 

380  PRINT  'THE  (LOLER!  FRINGE-SLOPE  UARIAELE  50(1)  IS  (  dl/dP.  )/ABS(  d2/dR  1' 

393  ■ - 

420  ■ - and  IliFUT  DATA  AT  EACH  INTERUAL - 

410  PRINT  'THE  INTERUAL  INDEX  I.  RHOIAL  (LEFT)  BOUNDARY  POSITION  RHI),(NEG.) 

AXIAL  POSITION  21  (  I  )  ,FSI,MGE-3L0PE  UARIABLE  33(1)  AND  THE  LOWER.  CENTRAL  AND 
UPPER  fringe-pair  SPACINGS  D 1 ( I  )  .02( I > ,03; I  )  APE  READ  IN  FOR  EACH  DATA  LINE.' 

420  . 21(1  )  IS  the  left-side  AXIAL  POSITION  OF  THE  LOWER  SHIFTED  FRINGE. FOR 

EACH  J  -  10  intersection  (WITH  SOdl-l).  21  (  I  )--D2(  I  1 ,  RATHER  THAN  . . 

420  FOR  I-:  TO  13 

440  READ  I  .R1  (  I  )  .21  (  I  I  .S0(  I  )  ,01  (  :  )  .02(  I  )  .D3(  I  i 
450  NEXT  l' 

460  input  'AT  THE  OUTER  RADIUS  C-IOH)  R  1  ,2  1  ,33  ,0 1  .02  ,ANO  03  ARE  '  i  R  I  (  10+ 1  )  , 

21 ( lO+l  ;  ,S0; I0+I  )  .01 ( I3  +  I  )  ,D2( lo+l ),D2; I3‘I  ) 

473  - - - - - - - - - - 

463  ■ - FIND  AXIAL  FRINGE-ST  POGIT'O-’IS  AT  I.'JTE.RU.-TL  BDUNCARIES - 

423  FOR  I  -  I  TO  10- 1 

503  G.  I  1  -  (C3(  :  )  -  0:  '  I  )  )/liO 

5  ;  .2  HI  ■  )  .  (  02  '  I  )  -  0  ;  .  I  '  >/  2C 
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s:ei  '  -  (oi ( i >  +  o:! n v:® 

533  8<  r  >  -  *  q.HiDi/fi 

543  Cl  I)  -  (3«u(1 1  -  .H(  T  ;  )/s 

550  NEXT  I 

560  POR  r-l  TO  I0+I 

570  XI  I .0  >-0 

530  Ztr,0)*ZH;) 

590  NEXT  I 

600  FOR  I  -  1  TO  10>l 

610  FOR  J  -  1  TO  9 

633  yil.J)  -  8(1)  4  (2*J  -  l)»CtI)/3 
630  XII.J)  •  J»(fl(I)  +  <  J-;  )«y(  I  ,J  )/3  ) 

640  Z( I ',; )  «  zi  1 1)  +  x< I ,j ) 

650  NEXT  J 

663  2(  1,10)  -  ZI  ( I )  *■  OZl : ) 

670  NEXT  I 

680  ■ - 

590  ■ - RIGHT-H/'NO  INTERUAL  BOUMDiiRY  VALUES - 

700  FOR  r  *  I  TO  10 

710  FOR  J  *  I  TO  10 

720  IF  Z'd*!)  •  0  OR  ZKI*-!)  -  -DZCtl!  THEN  750  ELSE  730 
733  XKI.Jl  -  X<I  +  1,J) 

740  SOTQ’tSO 

753  IF  ZI(I‘I>  -  0  then  763  ELSE  780 
763  XKI.J)  >  (01  (  1*1  )/C;(  I  )  )*X(  I  .J  ? 

773  SOTO  7S3 

783  XKI.J!  -  (03(  If  1  i/CZ' I  1  )*X(  I  ,J  ) 

790  IF  ZKIfl)  •  0  THEN  300  ELSE  830 
830  Z6(  I  .J  •  XI  (  I  ,J  )  -  01(1*1) 

810  GOTO  950 

830  IF  ZKI»!!  •  -03' 1*1  ;  THEN  930  E-SE  953 
833  ZS>  I  .J  )  •  XKI.J) 

943  SO'O  960 

853  ZS(  I  .J  )  -  ZI  (  I  +  l  )  f  XKI.J) 

860  NEXT  J 
873  NEXT  I 


880  . - - - - - - - 

893  ' . F;nO  hits  (FRI.NC-ELET  INTERSECTIONS  WITH  OR  STRADDLES  OF  CUT) 

903  . --RADIAL  POSITION  IN  INTERVAL  IS - 


910  FOR  I  -  I  TO  10 
920  FOR  J-  1  TO  9 

933  Y(I.J)  -  -Z(  I  ,J  )*(RI  (  Ifl  )  -  R1  (  I  ) )/(  ZG(  I  ,J  )  -  Z(I,J)) 

940  NEXT  J 
953  NEXT  I 
960  H  -  I 

970  FOR  I  •  I  TO  10 

980  IF  ZUI)  -  0  THEN  990  ELSE  1010 

990  J  -  0 

1000  SOTO  1183 

1010  IF  ZKI)  -  -02(1)  THEN  1023  ELSE  1043 
1020  J  -  10 
1030  GOTO  1 180 
1040  J  ■  I 

1053  IF  ABS(ZS(I,I)  -  Z(I.I))x.00l  THEM  1070  — FRINGELET  I  PARALLEL  TO  CUT-- 
1063  IF  ZS( I  .1  )>Z( I  ,1  )  ANO  ZS(I.I)>0  THEN  1180 
1370  FOR  J  -  2  TO  9 

1080  IF  ABSlZSd.J)  -  Z(I.J))<.001  THEN  1110  -FRINGELET  .(  PARALLEL  TO  CUT- 
1090  IF  J-3  goto' 1110 

1100  IF  Yd.J)>0  ANO  Y(  I  .J  XRl  (  I*l  )-RI  (  I  )  THEM  1130  '—INTERSECTION  TEST-- 
1110  IF  Z(  r.J  )/Z(  I  ,J-1  )<3  ANO  ZG(  I  .J  )/Z6d  .J-l  )<0  THEN  1263  '—STRADDLE  TE3T-- 
1  120  NEXT  J  '  ' 

1130  NEXT  I 
I  143  goto  1323 

I  1C3  - - - 

IIRX  - COUNT  AND  LA9E'_  HITS - 
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r.''3  ■ - INTERSECTION  COUNTER  ECULOUS - 

I  180  M  -  H  +  I 
11 S0  E^H)  -  0 
1200  I H  H  )  »  I 
IZIO  JKH)  -  J 

1223  IF  J-0  THEN  I01O  ELSE  1230 
1230  IF  J  -  10  THEN  l0iO  ELSE  12^3 

12- 10  IF  J-!  THEN  1073  ELSE  1123 

1250  ■ - STRADDLE  COUNTER  FOLLCLS - ^ - 

1250  H  -  H  +  I 
1273  E^Hl  -  .5 
1283  IHH'i  =  I 
1 2S0  J I <  H  )  -  J 
I  333  SOTO  I  123 

1310  - outer  inte.=s£;t:cn - 

1323  I  -  13 

1333  IF  zi;:3!  then  i3-ic  else  1353 

13- 13  J  >  13 
1353  -SCT:  13-3 
1353  3-3 
IS'O  H  -  H  -  I 

I  363  El  H )  -  3 
I  3S3  I ; ;  H  ;  >  I 
1-133  Jl'Hl  >  J 

I-1'3  P=:NT  -the  -‘1U7S£=  H  CF  HI'S  IS'.H 

1-123  2.7  F  H  ;  H  )  .R' H  1  ,RRf  H  .  ,3,  H  '  ,S  :  (H  ;  .321  H  1  ■H),'J2-.H1 

1-123  - CENTER  <R-31  HI* - 

1-1-13  r;  ;  )  -  0 
1-153  I!  <  1  I  -  1 

USO  E'  I  )  •  3 

l-l"3  IF  21(  1  1  >3  THEN  1-133  E-SE  1530 
1-163  J1  (  1  )  -  3 
USO  SOTO  15-13 
1533  FOR  J  •  1  TC  9 

1  =  13  IF  Z(  1  .2  1  -0  THEN  1523  ELSE  lE-13 
1  =  23  JK  n  -  J 
1530  30TC  16-10 

15- 13  NEXT  J 
I5E3  E'.  n  -  .5 

1553  IF  Z(I,I)  >0  THEN  1573  E-SE  1553 

1=70  J1  (  I  )  -  1 

1533  eOTO  1643 

15=3  FOR  J  -  2  TO  10 

16- 33  IF  Z(  I  ,J  )/Z(  I  .2-1  )<3  THEN  1610  ELSE  1630 
1613  J I ( I )  -  J 

1620  3073  1540 
1633  NEXT  J 

1543  - 

1653  ' - FIND  R.AOII  OF  FRINGELET  CRCS5INSS - 

16E3  FOR  I  -  2  TO  H  -  I 

1670  IF  ZlIKIl.JUIl)  <>  0  THEN  1653  ELSE  1730 
1630  IF  Eill  -  0  THEN  16=0  ELSE  1713 

16  =  3  R<  1)  -  RllIKD)  -  Z(  1 1  <  I  !  .J1  (  I  )  )»<R1  (  IK  :  !  ♦  1)  -  R  K  I  K  I  1  1  )/(  Z6v  I  K 
<  I  )  )  '  Z<  IK  I  1  ,JK  I  !  1) 

1733  SOTO  1743 

1713  Rill  -  .5«<RK;kI))  t-  RKIKtl  ^  I)) 

1723  5070  1740 
1733  ?<  I  >  -  RK  IK  I  11 
l~43  NEXT  I 
I  753  R-'  H  I  -  RK  IC-* ;  ) 

1753  'OR  r  -  I  rc  H 

I  T-.o  i.-'R  ;n:  -h  I  t :  a:  h  e  i  i  .r  -  ■ ;  i  e  ; . : : ■ ' :  > ;  k  -  ;  i  nc :  i 
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1780  NEXT  I 
1790  ’STOP 
1800  ■ - 


1810  ■ - INITIAL  FRINGE  CALCULATIONS - 

1820  - 

970  ' - SPAN-TO-INTERUAL  INDEXING - 

1830  PRINT  "THE  SPAN  BOUNDARIES  ARE  LOCATED  AT  I  •  I ,10  +  I .  AND  AT' 

18iO  ■ - 1  >  K  -1  FOR  THE  FOLLOUING  LOOP - 

1850  FOR  I  •  1  TO  KO-1 
1860  INPUT  'IZIK-l  )  -• i I3( I ) 

1870  NEXT  I 
1830  K  »  0 
1890  U2  -  1 

1900  ■ - 1  -  interval  for  FQLLCUING  LOOP - 

1910  FOR  I  -1  TO  I3(K0  -  1  ) 

1920  IF  t<r3<U2)  then  1940 


1930  U2  -  U2  t  1 
1940  K( I )  -  U2 
19S0  NEXT  I 

1960  FOR  I  =■  I3(K3  -  1  )  H  TC  10 
1970  K< I  )  -  KO 
1980  NEXT  I 

1990  ■ - 

2000  PRINT  ‘THE  LEFT-HAND  SPAN  BOUNDARIES  ARE  AT  RADII  '  ' - R2(  n  =  0 - 

2010  FOR  I  •  1  TO  XO 
2020  INPUT  -RS'  X  :  ;R:(  I  > 

2030  NEXT  I 

2040  RC'KOt:  1  -  R!(  lOM  ) 


2060  PRINT  ‘ThE  INITIAL  L3'-£R  FRINGE  NUNSEPS  ARE- 

2070  FOR  I  -  I  TG  KO 

2030  INPUT  -N!  (K  )  -•  lA.i;  I  1 

2090  NEXT  I 

2100  . — - - - - . - 

2110  ■ - - FIND  l.K.N.MI  at  hits - 

2120  ■ - - — - - - - . 

2130  PRINT  -THE  LEFT-SIDE  AXIAL  POSITIONS  <NE3.  )  OF  THE  LOUER  INITIAL  FRINGES 
ARE" 

2140  for  I  -  1  TO  KO 
21S0  INPUT  -ZC'.K  )  -■  iZ2(  I  ) 

2160  NEXT  I 

2170  • . — - - 

2180  PRINT  -THE  RIGHT-SICE  AXIAL  POSITIONS  (NE3.)  OF  THE  LONER  INITIAL  FRINGES 
ARE' 

2190  FOR  I-I  TO  KO 
2200  INPUT  ■Z4(K !-■ iZ4( I ) 

2210  NEXT  I 

2220  - 

2230  PRINT  -THE  LEFT-SIDE  AXIAL  POSITIONS  (POS.)  Z3(K)  OF  THE  UPPER  INITIAL 
FRINGES  ARE- 
2240  for  I  -  I  TO  KO 
2290  INPUT  -ZGIK  ;23(  I  ) 

2260  NEXT  I 

2230  PRINT  -THE  RIGHT-SIDE  AXIAL  POSITIDNS  (PCS.)  OF  THE  UPPER  INITIAL  FRINGES 
ARE- 

2290  FOR  I-  I  TO  KO 
2300  INPUT  -ZSIK  )-- iZ5(  I  ) 

2310  NEXT  I 


- initial  FRINGE  FRACTION  AT  EACH  HIT-- - 

2330  ' - - - UK  (O'.  AND  U2(>0)  ARE  AXIAL  POSITIONS  OF  LOUER  AND  UPPER 

FRINGES  AT  A  HIT,  WHILE  SI  AND  S2  ARE  THE  FRINGE  SLOPES  — . . 


2  39,1  Fill  .  -z;{i)/(z:(i)  -  zz'i', ) 
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23S0  FOR  I  -  2  TO  H-1 

2370  SKI)  -  (Z4(K(IKI)))  '  Z:(K(Il(I))))/(R2<K(tKI))  +1)  -  RZCKdKt)))) 

2360  S2(I)  -  (ZS(K(IKI)))  -  Z3(  K(  I K  I ) ) )  )/(R2<  K<  1 1  <  I ) )  +1)  -  RZiKCIKI)))) 

23SO  UHI)  -  ZitKCIKI)))  -  SKI)*(R2(K(I1(I1)  +1)  -  R(  I ) ) 

2400  U2(  I  )  -  Z5(K(IKI)))  -  S2(I)»<R2(K(II(I))  +1)  -  R(  I ) ) 

2410  F<n  -  -UKI)/<U2(n  -  UKD) 

2420  NEXT  I 

2430  F(H  )  -  ( 1  -  S0( 10  )  )/2 
2440  FOR  I  -I  TO  H-1 

2450  PRINT  ‘H.Sl  ,S2,U1  ,U2>'iI:SKI)iS2(I):UUI)iU2<I) 

2460  NEXT  I 
2470  'STOP 

2480  - 

2430  ■ - FINO  the  phase  SHIFT  P(  I  )  AT  EACH  HIT - 

2500  FOR  I  -  I  TO  H-1 

2510  PlI)  -  2»(N(L(I1(I  )))  +  (JKI)  -  E(I))/10)  -  2»(  MK  K1  1 1  ( I ) ) )  +  F(  I ) ) 
2520  NEXT  I 
2530  P(H)  -  0 
2540  'STOP 

2550  • - 

2560  ■ - SORT  AND  INDEX  FOR  INCREASING  RADIUS - 

2570  FOR  I  -  1  TO  H 
2530  s; I  )  -  I 
2530  RR( I  )  -  R( I ) 

2630  NEXT  I 

2610  NR  -  H 

2623  IS  -  3 

2630  NR  -  NR  -  1 

2640  FOR  I  •  I  TO  nR 

2650  IF  R(I)  <  Rd  +  l)  GOTO  2733 

2663  IS  -  1 

2670  RH  •  R(  IM  ) 

2633  R< I+:  )  -  RC I  ) 


2633  R( 1)  -  RH 
2703  SH  -  S(  I  +  l  ) 

2710  s( IK)  -  sa > 

2720  S( I  )  -  SH 
2730  NEXT  I 

2740  IF  IS  -  1  GOTO  2620 

2750  ■- . - . - 

2763  ■ - 

2770 
2783 
2790 
2800 
2810 
2820 


-final  results- 


iS(  I  )  .N(L(  IKS(  I  ) )  )  )iril(K<  II  (S(  I  )  )  )  )  .F(S(  I  ) ) 


FOR  I  -  1  TO  H 
LPRINT  -H  N.fll  F> 

NEXT  I 
’STOP 

FOR  I  -  1  TO  H 

LPRINT  -H  ,E  .1  ,J  ,R  ,P‘  :S(  I  )iE(S(  I  )  );IUS(  I  ))i  JKS(  I  )  )  .RR(S(  I  )  )  ,P(S(  I  )  ) 


2830  NEXT  I 
2840  STOP 

2853  ■ - 

2863  ■ - 


-OATA- 


2973 

OATA 

1  ,0  19  .-I  .  .4 

..4  ..35 

2880 

DATA 

2  , .87  . -.22,-1  . 

.38  . .38, .38 

2890 

OATA 

3.1  .78, -.3, -I  , 

.4 . .4  ,  .4 

2900 

OATA 

4, 2. 52. -.26.-1 

..38  , .33, .32 

2910 

OATA 

5, 3. 63, -.23,-1 

..36, .38, .34 

2920 

OATA 

6. 5. 02. -.25.-1 

..36  ..35, .30 

2930 

OATA 

7,  6. 50. -0.0, 

-1  ,.37, .32.. 

2940 

OATA 

8  .  7.10  .-.08 , 

-1  ,.38. .30.. 

2950 

DATA 

9,  8. 40. -.14. 

-1  ,.3a,.30,. 

2960 

OATA 

10  .  9.02  .-0.0  . 

-1  ,.28. .30  .. 

2973 

OATA 

11,  9.28 .-  .32  , 

-1  ,.30, .32.. 

2990 

OATA 

12.  9. 45. -.07. 

-1  ..31  ,.33,. 

2990 

OATA 

13,  9. 57, -.13, 

0,.27,.33.. 

3030 

DATA 

14  ,  9.90 ,-. 12 

1  ..24, .30,. 

30:0 

OATA 

IS  10. 18  -.25 

1  . .23 . .25  ,  . 

20 


3020  DATA  16,10.37.-. 14.  I  ..20  .  .24  ,.33 
3030  DATA  17, 10. 69. -.10,  1  .  .29  .  .22  .  .32 
3040  DATA  18, 10. 92, -.10,  0  ..30  ..22  ..33 
3060  DATA  19,11  .61  .-.13,-1  ,.32.. 27,. 29 
3060  DATA  20, 12. 62. -.17. -I  ,.30. .28. .29 
3070  DATA  21 .13. 68, -.20,-1  ,.29. .29. .30 
4000  END 


Final  results/  giving  hit  number  (H),  straddle  index(E)/  interval 
index  (I)/  fringelet  index  (J),  and  radius  (R)  and  phase  (P)  at 
each  hit. 


H.E.I  .J.R.F 

1 

.5 

1  5 

0 

3.7 

H  ,E  ,I  .J  .P  .P 

T 

0 

1 

5 

.2G9G1G7 

3.678824 

H.E.I  .J.R.F 

3 

0 

2 

G 

.9770586 

3.5G0872 

H  .E  .1  .J  .R  ,P 

4 

0 

2 

7 

1 .504242 

3.523936 

H.E.I, J.R.F 

5 

0 

3 

7 

2.253S 

3.195322 

H.E.I  .J  .R.F 

G 

0 

4 

G 

3.G7S732 

2.396723 

H.E.I  .J.R.F 

7 

0 

5 

G 

3.70243 

2.3854! 1 

H.E.I  .J.R.F 

8 

0 

5 

7 

5.008705 

2.031905 

H  .E  .1  .J  .R  .P 

9 

0 

G 

8 

5.504912 

2.021648 

H.E.I  .J  .R.F 

10 

0 

B  9 

G .000786 

2.01 1532 

H  .E  .1  ,J  .R  .P 

1  1 

0 

7  0 

6.5 

2 

H  ,E  .1  .J  .R  .P 

1  2 

0 

7  2 

G . 997422 

2.  174922 

H.E.I  .J.R.F 

13 

0 

8  3 

7.4S274S 

2. 1G43G3 

H  .E  .1  .J  .R  .P 

14 

0 

8  4 

8.  1235G3 

2.GE535G 

H  .E  .  I  .J  .R  .F 

15 

0 

9  5 

S.4G382I 

2. 106678 

H  .E  .1  .J  .R  ,P 

1  G 

0 

9  G 

8.53581 

2 . 256 1 95 

H.E.I  .J  ,R  ,P 

17 

_  0 

9  1 

8.697731 

2 . 40555 1 

H.E.I  .J.R.F 

18 

£ 

5  6 

8.80672B 

2.5=6233 

H.E.I  .J.R.F 

1  9 

0 

9  9 

8.913634 

2.707768 

H  .E  .1  .J  ,R  ,F 

23 

0 

10  0 

9.020Q0I 

2.859723 

H.E.I  .J.R.F 

21 

0 

1  1  2 

9.429331 

3.074239 

H.E.I  , J.R.F 

22 

0 

12  3 

9.543399 

3.222897 

H.E.I  .J.R.F 

23 

0 

12  4 

9.657842 

3.371  1  12 

H.E.I  .J.R.F 

24 

c 

D 

13  5 

9.785 

3.413574 

H  .E  .1  .J  .R  .P 

23 

0 

1 4  4 

9.919374 

3.252772 

H  .E  .1  .J  ,R  .P 

27 

0 
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VI.  APPENDIX 


Interpolation  of  the  non-unif ormly-spaced  shifted  fringes  is 
accomplished  by  dividing  the  axial  space  between  the  two  fringes 
which  bound  the  cut,  within  an  interval,  into  ten  non-uniform 
spaces  by  introducing  decimal  "f ringelets" .  The  lower  bounding 
fringe  is  taken  as  the  J-0’  fringelet  and  the  upper  bounding 
fringe  is  taken  as  the  3=10  fringelet.  Refer  to  Fig.  2.  Let 
Dl ( I ) ,  d2 ( I ) ,  and  D3(I)  denote,  respectively,  the  axial  fringe- 
pair  spacing  at  the  (left)  interval  boundary  of  the  lower, 
central  (bounding  the  cut),  and  upper  fringe  pairs.  It  is 
necessary  to  determine  the  distance  D(I,J),  for  J=1  to  10, 
between  the  J-1  and  J  fringelets.  First  require  that  the  sum  of 
the  fringelet  spacings  be  the  central  fringe-pair  spacing. 


I , J)=D2( I ) 
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Information  about  the  non-uniform  adjacent  fringe-pair  spacings 
is  included  by  requiring  that  the  lower  fringelet  spacing  D(I,1) 
be  one-tenth  of  the  average  fringe-pair  spacing  for  the  lov/er  and 
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central  fringe-pairs  and  that  the  upper  fringelet  spacing  D(I,lf() 
be  one-tenth  of  the  average  fringe-pair  spacing  for  the  upper  and 
central  fringe  pairs. 


D(I,1)-(D1(I)+D2(I)  )/2jef  A2 

D{I,1J0)-(D2(I)+D3(I)  )/2j0  A3 

In  order  to  meet  these  three  requirements,  D(I,J)  must  be  at 
least  quadratic  in  J.  It  is  assumed  that 

D( I , J)=A( I )+B( I ) ( J-1 )+C( I ) ( J-1 ) 2 .  A4 

The  coefficients  A,  B,  and  C  are  determined  by  using  Eq .  (A4)  in 
Eqs.  {A1-A3).  It  is  immediately  seen  that  A(I)  is  D(I,1)  and  is 
thus  given  by  Eq.  (a2).  The  other  coefficients  are  evaluated, 
after  a  little  algebra,  as 

B{ I )=(-19*G( I )+9*H( I ) )/8  A5 

C(I)=(3*G(I)-H(I))/8  A6 
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where  G{I)  and  H{I)  are  defined  by 


G(I)»(D3(I)-D1(I)  )/18j?  A7 

H(I)*(D2(I)-D1(I)  )/3if  •  A8 

In  calculating  the  intersections  of  the  fringelets  with  the 
cut,  it  is  necessary  to  know  the  distance  X(I,J)  between  the 
lower  bounding  fringe  and  the  J  fringelet.  This  distance  is  the 
sum  from  J'=l  to  J  of  D(I,J').  Evaluation  of  the  sum  shows  that 

X( I , J)=J*{A( I )+( J-1)*V( I, J)/2)  A9 

where,  V( I , J ) aB ( I ) + ( 2J-1 ) *C( I )/3  •  AlO 
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